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1. Introduction 

The dynamics of relaxing quantum systems in sufficiently strong laser fields, which are 
neither stationary nor monochromatic, is an important aspect of the theory of interaction 
between atoms and photon fields. This problem is of particular practical significance 
because of dramatic advances in the precision of measurements presently attainable 
in spectroscopic experimental studies of hydrogenic and few-particle atoms. Among 
the simple atoms that attracted attention in the last few years one should mention: 
hydrogen and its isotopes [1, 2, 3, 4, 5], positronium [6, 7], denoted (e"*" — e~), muonium 
(see [8, 9, 7, 10, 11] and references therein), denoted (/i"*"— e~), and helium atom [12]. The 
measurements of the ground state HFS and the Lamb shift are of particular interest for 
these systems. High precision spectroscopy of low-lying transitions in hydrogenic atoms 
offers a unique opportunity to test QED calculations and to refine both the fundamental 
constants and various properties of respective nuclei. This urges further theoretical 
developments intended to describe the well established photon-induced processes in 
simple atoms with much higher precision than techniques used so far are able to provide, 
in order to allow a proper interpretation of newly available experimental data. In the 
present paper the stepwise 3-photon ionization of muonium in a non-monochromatic 
laser field with a time-dependent amplitude is studied. 

It has been recognized ah^cady a long time ago (sec [13] for a more detailed 
discussion) that the IS — 2S transition offers unique opportunities for high precision 
spectroscopy due to the narrow natural line width, r^2j'^\ of the 25'-state. In the 
hydrogen atom, for example, F^^^*^ 1.3 Hz which enables the quality factor of 
~ 10~^^ to be achieved already today, as well as suggests future resolutions of order 
of 1 part in 10^^. Experimentally, the l^S" — 25" transition can be induced Doppler-free by 
absorbing two photons from two counter-propagating laser beams. In particular, this 
scheme is presently used at the Rutherford Applcton Laboratory [8] for determination 
of the IS — 2S energy separation in muonium to the 1 MHz accuracy. This study offers 
an opportunity to improve the present knowledge of the muon/electron mass ratio and, 
thereby, of the muon mass itself. In the course of this investigation, two counter- 
propagating pulsed laser beams with almost identical intensities, I{x,y,t), /max ~ 
10^ W/cm^, and the same photon energies, hi/L ~ 3/16 a.u. (Al = 244 nm) are used. 
Thus, the frequency i^l is such that it is tuned into resonance with the 2-photon Doppler- 
free IS 25" transition. Unlike the hydrogen case however, the IS — 2S transition 
in muonium is hard to be observed directly, i.e. by detecting radiation emitted as a 
result of the 2iS'-statc dccxcitation. This is due to the fact that the appropriate line 
intensities in muonium happen to be weaker by several orders of magnitude, owing to 
much lower densities at which muonium atoms can be produced. Furthermore, the 
latter reason necessitates the use of higher laser intensities (by a factor of 10^, at least) 
that are required for a reasonable signal strength. This demands a pulsed rather than 
a CW laser source (as in the case of hydrogen [5]) to be employed. The 25'-state of 
muonium is detected therefore via its photoionization by the third photon absorbed 
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Figure 1. (a) - The scheme of the new 1S-2S experiment in muonium. (b) - The set 
of states to model the stepwise 3-photon ionization. 



from one of the laser beams [8] . Schematically, the experiment in muonium is presented 
in figure 1(a) where hwi, 3/16 a.u. (i^L = oj-l/It^ ^ 10^ MHz) stands for the resonant 
energy (frequency) which drives the 2-photon 15" — 25" transition. 

Despite the fact that the laser frequency can be calibrated, in principle, to a 
rather high accuracy, there are a number of systematic error sources each of which 
proves to be essential for precise determination of the 15* — 25* energy separation in 
muonium. In this paper we report the results of the theoretical study intended to allow 
for one of the most important among these effects, that is, the time-dependent frequency 
variation (chirp) of the laser field. This phenomenon arises due to the rapid refraction 
index variation of the laser media and is hardly avoidable, especially with powerful 
pulsed lasers, unless some technical developments are made in order to compensate the 
frequency alteration [8, 12]. As will be demonstrated below, the chirped laser signal 
immediately makes its appearance in appreciable spurious shifting and broadening of 
spectral/photoionization lines, as well as leads to a noticeable distortion of the line 
shapes. To estimate the scale of the effect considered, one should note, for example, 
that the chirp-induced shift of the center of the ionization line is roughly equal to a 
characteristic magnitude of the chirp itself. For the 1S — 2S experiment in muonium, the 
latter can typically amount to 5/27r ~ 10 . . . 100 MHz, thus leading to the relative shift 
of the line's center by S/coi, ~ 10~^ . . . 10~^. The effects of such an order are usually 
completely neglected in atomic physics. Nonetheless, these happen to be crucial for 
studies where an absolute accuracy of 1 MHz is anticipated. It is therefore of primary 
practical importance to work out a relatively simple theoretical scheme which could allow 
one to describe certain types of multi-photon ionization phenomena in hydrogenic atoms, 
where an arbitrary time variation of the frequency of the laser field is accompanied by 
an arbitrary time modulation of its amplitude. 

Although the problem of the hydrogenic systems' interaction with monochromatic 
resonant laser field with/without the amplitude modulation is investigated quite fully by 
now (see [14, 15, 16, 13, 17, 4] and references therein), there have been much fewer studies 
where both the amplitude and the frequency of the laser signal vary with time. Even 
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though, the authors confine themselves in most cases to either 1-photon (rather than 
2-photon) resonant transition [18], or treat only somewhat special forms ("abrupt step", 
"exponential field pulses" etc.) of the field modulation [19]. Under these circumstances, 
it seems very desirable to reconsider the problem as a whole, while adapting it specifically 
to the conditions of the IS — 2S experiment in muonium. 

2. Theory 

2.1. Basic equations 

We consider, without a loss of generality, a set of three levels: |a), \b) and |c), of a 
hydrogenic atom with the charge Z of its nucleus, such that their one-particle energies 
satisfy the inequality: £„<£;,< ^c, this rather general setup is shown in figure 1(b). It is 
supposed that the given 3-level system is exposed to two counter-propagating (along the 
;z-axis) laser waves with equal time-dependent (circular) frequencies, uj{t) = uji^ + j(f){t), 
polarizations vectors, e, and wave vectors, ki(t) = —k2{t) = k(t), \k(t)\ = auj(t), where 
a = e^/ (he) ^ 1/137 is the fine structure constant and c denotes the speed of light. The 
stationary part, uji^, of the field frequency a;(t) is assumed to be such that it is in the 2- 
and 1-photon resonances with the pairs of states, (|a), \b)) and (|6), |c)), respectively. 
This implies that 2^1/]^ = cu]^ — ujb^a/'^ — ^c,b with cuj^i = Sj — ei, — a,b,c) being 
the differences of the one-particle energies. In addition, it is supposed that the 
transition is forbidden between \a) and \b) and is allowed between \b) and |c). 

Under the resonance conditions assumed above, 3-photon ionization of the system, 
that is, I a) |c) transition, occurs predominantly as a 2-step process: (i) a 2-photon 
resonant absorption from the state \a) into the state \b), followed by (ii) a 1-photon 
resonant transition between the levels \b) and |c). This simple physical picture is a 
consequence of the evident estimate of the probability for an atom to absorb 3 photons 
of equal energies (the atomic units, e — H — nie = I, are used throughout the paper; in 
this system, the speed of fight is equal to c ^ 137): 

W^i2(res) oc ^<j(res)<i)(res). (1) 

Here, wl"^^ (res), W^^,j{res) are the 2- and 1-photon resonant transition probabilities. 
Fa is the total width of the state | b) , and Acu denotes some characteristic mean energy 
which a two-fold summation over intermediate states in the exact expression to define 
Wi.a' (res) (see [20] for a more detailed discussion) can be reduced to. Equation (1) readily 
follows then from the estimate, Ta/ Acu <ti 1, which is fulfilled for most states of atoms 
in moderately strong laser fields, except extraordinarily short-living ones. A situation 
when I a) stands for the ground-, \b) for the 25*-, and |c) for the continuum eP-states, 
respectively, is of our primary concern here. In this case Auj^ ~ 1 and T^^°'^/aJ ~ 10"^ 
since Fg*"*"* = F^J"^^ + r2^'^°*'* ~ 10~^ is dominated by the photoionization rate of the 2S- 
level {T^2^"^^ — 10~^ ^ ^^Z^^^ — 10~^^). Accordingly, we will assume that the natural 
line widths of the \a)- and |6)-levels are both equal to zero: ri"^^*^ = p(nat) _ 
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straightforward generalization of the latter condition, that might be required for the 
treatment of excited states, can be achieved by adding appropriate imaginary parts to 
the energies of the levels [21]: Sa ^ £a — iTa^^^^ /2, Sb ^ £b — ir[°'^*-'/2- 

In the semiclassical approximation, which happens to be accurate enough for our 
purposes, the electric fields inducing the 2- and 1-photon transitions can be taken in 
the form 

Ei,2{r, t) = ^eEi,2(t)t/i,2(r") exp {i (fci,2(t)r - uj{t)t)} + c.c. (2) 

Here, the real functions Ui{r) = U2{r) = U{r) and Ei(t) = E2{t) = E{t) describe the 
(macroscopic) spatial inhomogeneity of the laser field and its time-dependent amplitude. 
These are related direct to the laser intensity I{r,t) which, along with the chirped 
(circular) laser frequency uj(t), is obtainable direct from measurements. For the 15* — 25* 
experiment with muonium, I{r,t) and Lv{t) can be well approximated as [22]: 

Hr,t) . ±W(r) . p^^^cp {-^ - ^ , (3) 

1 • 1 

uj{t) =ujL + -Ht)^^h + (f>{t)=uj'L + -S{l + erf{t/T)), (4) 

where — 2 ... 5 mJ is the energy within one laser pulse; fJt = 40 ... 65 ns, 

(Tr = 0.5 . . . 1.5 mm denote the temporal and spatial dispersions of the external laser field 
(in the XOY plane perpendicular to the direction of the beams' propagation). The full 
temporal width of the laser pulse is expressed then as r = 2^2 ln2 ■ = 100 . . . 150 ns. 
Parameter 6 = cjl(+oo) — u;l(— cxd) ~ 27r- (10 . . . 100) MHz defines the magnitude of the 
chirp 0(t), being actually detected in the experiment, and erf(. . .) stands for the error 
function [23]. 

As usually happens to be true in most practical situations, both the amplitude 
E{t) and the additonal phase (f){t) vary much slower with time than E{t) and u;{t){see 
equations (3)- (4)). In addition, the magnitude of the chirp is normally much smaller 
than the "base" frequency cuj^. It is therefore permissible to assume that u!{t) and E{t) 
are subject to the following general conditions [24, 25]: 

1 



\tuj{t) \ -C \uj{t) -culI 



|0WI«^L (5) 



E{t) <^\u;j^E{t)\. (6) 

Within the framework of the conventional time-dependent perturbation theory [21], the 
exact time-dependent wave function of the quasi-stationary state ipnit) in the field is 
sought in the form 

V'n(i) = Ec,,„(t)V'fe--'=*, (7) 

k 



where Ck^nit), being the amplitudes to find an atom in one of the unperturbed states 
'4'k'\ satisfy the following well known system of coupled differential equations: 

Ck,n{t) = -i I] Vk,s{t)cs,n{t)e-''''-'. (8) 
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Here, Vk^s{t) denotes the matrix element of V{t) being the operator of the "particle-laser 
field" interaction. Due to the given below reasons, V{t) can be taken in the form 

V{t) ^ -E{r, t)'d^ -^E{t)U{r)e-'''^*^*e - d + c.c. (9) 

where d is the dipole operator of an electron (or a muon) and E{r,t) denotes the 
electric field (2) seen by the atom in its rest frame. For particular laser parameters of 
interest, the accuracy of the above approximation follows from the evident estimate: 
|VC/(r)| ~ l/(7r 10~^ -C \k{t)-r\ ~ 27r/AL 10~^. This implies that the contribution 
of the quadrupole terms oc £^(i)(VC/ (r) •</), which originate from the spatial macroscopic 
inhomogeneity of E{r, t), is expected to be a factor 10"^ smaller than the contribution 
of the (dipole) terms in (9). Almost the same estimate holds also for corrections 
arising due to the quadrupole component oc E{t)U{r)\k{t)\{e • dY of the photon field 
itself. This is despite the fact that the states \a) and \b) are coupled, in accord with 
the initial assumptions of our model and equation (9), in only the second order of 
the perturbation theory in V{t), i.e. by taking into account at least the 2-photon 
absorption/emission. Indeed, unlike the second order dipole coupling, the 1-photon 
quadrupole \a) <-> \b) transition would apparently be off resonance. This would makes 
its appearance in the extra exponential factor ~ e^'^^* in the amphtude. If K,6 <S 
as it is assumed throughout the paper, then the terms of such a type are known to be 
almost negligible [26] . In our case, in particular, their relative contribution is of the order 
of 1/(c<jlt) ~ 10~^, as can be shown by averaging the transition amplitude over some 
time interval: —T < t < T c:^ r. The same "off-resonant" argument is fully applicable 
to magnetic transitions as well, thus elucidating why the Ml- and £^2-contributions 
are suppressed by a factor ~ 10^^ as compared to the dipole one. 

By virtue of (8)- (9), one can readily retrieve, after some conventional algebra, the 
following form for effective quasiclassical one-particle operators which couple two given 
states of reference [17, 13]: 

vS(r,t,'^L) = ic/^(r)E'(t)a-'^""^'^^")'^ -^I''><''l'^-^l''> , (10) 

Vi^l{r,t,u^) = i^^(r)E^(t)e^'(-^^+'^W) >f ^"'^ ' "^'^^^^'^ ' , (11) 

All summations are performed here over complete set of one-particle atomic orbitals. 
Equations (lO)-(ll) describe the real transitions, whereas (12) defines the AC Stark 
shift of the level \b) (or |a), if 6 — > a) which arises due to interaction of the atom with 
the incident radiation. Note that the sum in (12) coincides, by definition [21], with the 
dipole dynamic polarizability of the state at the frequency uj. Also, it should be 
pointed out that only the leading terms, as given by equations (10)-(12), were retained 
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in the course of derivation. The most important among neglected terms are of the form: 

2 



U\r) 



.E{t) 



-iO-i^)t-<l,(t)] 



+ c.c. 



m U\r). 



Apparently, this contribution is smaller than V^^j {r^t^uoi), = a,b) by a factor of 

E{t)/{uJi^E{t)) ^ ~ 1/(c^lt) ~ 10-^ as follows from (3)-(4), in accord with (5)-(6). 

As was assumed above, there are no other discrete resonant atomic levels (either 
1- or 2-photon) except |a), \h). This enables the system (8) to be reduced to these two 
states only §, while describing the 1st stage of the two-step 3-photon ionization process 
under consideration. Indeed, an inclusion of the off-resonant states would give rise in 
the right-hand side of equations (14)-(15) to the rapidly i-varying terms ~ Q^'^'i^*, After 
temporal averaging (see above), these happen to be of order of 1/{uji,t) ~ 10~^ , and 
can therefore be omitted to the accuracy adopted in this study. On introducing the new 
functions, 

1 



-a,a(^) 



-ia(t) 



C{t) = Cfo,a(t)exp<^ i 



UJ 



^b,a ]+a{t) + 24>{t) 



(13) 



and replacing the matrix elements in (8) by equations (10)-(12), one arrives at the 
following system of equations: 



47r 

a{t) ^—I{r,t) [Da,a + Da,bC{t)], 



J^b,b J 

c 



(r, t) - (attAu + a{t) + 20(t)) 



(14) 
(15) 



where 27rAz/ = c<jl— c^6,a/2 is the time-independent frequency detuning off the resonance, 
I{r,t) = cE'^ {t)U'^ (r) / Sn is the intensity of the single laser beam and Dj^i, {j,i) = a,b 
denote the effective matrix elements given by the sums in (10)-(12). It has been taken 
into account here that the 3-level system is actually probed by a superposition of 
two laser fields (2), rather than by a single one; this results eventually in doubling 
all matrix elements. In addition, we have neglected a contribution due to the 2-photon 
absorption from each of two laser beams. This contribution comes from the second-order 
interactions of an atom with only one of two light waves and makes its appearance finally 
in two Doppler-broadened terms, i.e. depending on a f ^-component of atomic velocity in 
the laboratory frame, in the expression to define ionization probabilities. However, after 
averaging over v^, as would be required within the framework of a systematic approach, 
these terms turn out to be small compared, at the center of the photoionization line 
[27, 28], with the Doppler-free ones which originate from (14)-(15). These arguments 
provide the ground for the approximation made, although the effects caused by atomic 
movement in the media are beyond the scope of our present consideration. 

It is relevant to point out here that the matrix elements, Dji, = a,b, in (14)- 

(15) are c<;l ^dependent. In particular, Dbb has a non-zero imaginary part. The latter 

§ This implies that the Coulomb degeneracy of the levels with the same principal quantum number, 
say \ns) and \nd), is ignored, since we are mainly concerned here with the lowest IS- and 2S'-states. 
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C{t) (16a) 



allows for a non-zero photoionization rate of the |6)-statc, due to the action of the 
photon field, since cul is supposed to exceed the photoionization threshold /{,. Indeed, 
as was noted above, D^ h is proportional to the dynamic (tensor) polarizability ^^■'(cjl) 
of the level \b), which is complex-valued if cjl > h (see our recent paper [31] for a more 
detailed discussion). ^a^J^uj) and Qa'i^{Lj) define then the AC Stark shift of the level (see 
equation (12)) and its decay probability (oc iHcP) via the single photoionization. This 
result is in agreement with the optical theorem [20]: al"'\LJ-i,) = Anaooi^'^a,^ {ui,) , where a 
is the fine structure constant and a^^ (cui,) denotes the total photoionization cross section 
of the state \b). Numerical calculations show that the single photoionization happens to 
be the main channel of the |6)-level depopulation. It is therefore crucial for the correct 
description of the time-evolution of the 3-level system that this depopulation mechanism 
is taken in fact into account in (14)-(15) through the matrix elements Dj^i, = a,b. 

An allowance for the 2nd step of the process considered can be made by augmenting 
the system (14)-(15) with an additional equation that explicitly describes the \b) — \c) 
coupling through the single photoionization of the state \b). Given by (16c), this extra 
equation is nothing else but a differential form of the conservation law which determines 
the probabihty balance between the states of reference. Moreover, equations (14)-(15) 
can be decoupled by substituting a{t) of (14) into (15) to finally get 

ATI r 4:71 / • \ 

C{t) = -iDh,a—I{r, t) + i {Da,a ' D,,,) —I{r, t) + 2 (271 + 0(t) j 
+ iD,,Al{r,t)C'(t), 

An 

a(t) = —I{r, t) [Da,a + D,,aC{t)] , (166) 

Wc{t) = —at'\iJi^)I{r,t) \C{t)fexp{2Qa{t)}, (16c) 

where Wc{t) is the probability for an atom to be ionized by the time t, via absorption of 3 
photons of equal energy. Under the actual experimental conditions [8], initial population 
of the level \b), i.e. before the laser pulse is shot, is usually rather low. This is because 
the major fraction of the muonium atoms is produced in the ground state \a) where the 
population is proportional to |ca_a|^ = exp {2S>Q;(t)}. In view of this physical picture and 
by virtue of (13), it is natural to impose the following initial conditions on the functions 
a{t), C{t), W,{t): 

C(-oo) = 0, a{-oo) = 0, We(-oo) = 0. (17) 

These enable Wc{t — +00, Ai/, A^^, r; [0]) to be determined uniquely. Being of our 
primary concern, this quantity constitutes the required 3-photon resonant ionization 
probability of an atom at infinitely large positive times, i.e. when the interaction 
between the single laser pulse and the system has already ceased. As a function of 
the laser frequency detuning Az/, Wc(t = +00, Au, A^,r; [(f)]) describes, and will be 
therefore used as a synonym of, the 3-photon ionization line profile/shape. For brevity, 
the latter will be denoted as W^lAu, S, A^, r; [(f)]) where the argument [0] is introduced 
in order to indicate explicitly that the line profile depends on the form of the chirp. 
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In particular, if = 5 (1 + erf (i/r)) /2 then W^{Aiy, A^,r; [(f)]) will be denoted as 
simply W^{Au, A^, S, r). 

Equations (16a)-(16c) are well suited for numerical calculations with arbitrary 
functions defining the spatial and temporal distributions of the laser pulse. In addition, 
the given system is rather convenient for analytic treatment, as will be demonstrated 
in subsection 2.2. In particular, one can develop further perturbative expansion of 
equations (16a)-(166) in terms of the intensity I{r,t). This yields analytic formulae 
which, after setting ^(i) = 0, agree with those formerly derived in [17] within the 
framework of a 2-level model, for unchirped laser signals. Despite a relatively simple 
form of the appropriate equations, however, a not very straightforward numerical 
integration is still required in order to calculate, for example, photoionization line profiles 
in most practical situations, let alone the fact that the result of the work [17] is valid for 
rather weak laser intensities only. Note that equations (16a)-(16c) arc free from the latter 
restriction. Furthermore, unlike our approach employing the probability amplitudes, the 
problems similar to those considered in this paper are usually treated (e.g. see [17, 12]) 
in the formalism of the density matrix [29]. Under the physical conditions adopted here, 
these two approaches are equivalent. To facilitate the adequate comparison however, it 
is relevant to present here an alternative form of (16a)-(16c), by rewriting the system 
in terms of the density matrix elements: Pa,aif) = |ca,a(^)P, Phfiif) = \cb,a{t)\'^ and 
Pb,a{t) = c&,a(^)<„(t)exp{-i 2{2TTAut + (f>{t))}, denoted as pb,a{t) = -f{t) + ip{t). This 
yields: 

Stt Stt 
Pb,b{t) = —^Db,bl{r, t)pb,bit) - —Db,a I{r, t)(5{t) (18a) 

m = 2 (27rAi/ + 0(t)) 7(t) + ^/(r, t) {(3?D„,„ - m.b) 7(^) 

47r 

+ {^D,,a + $>A,6) m) + —Db,a Hr, t) {pb,b{t) - Pa,a{t)) W 
Pa,a{t) = ^^D,,a I{r, t)pa,a{t) + ^Db,a I{r, t)P{t) (18c) 

7(t) = -2 (27rAi/ + 0(t)) (3{t) + ^/(r, t) {(^,6 - ^Da,a) m 

+ (m,a + m,6)7(^)} (18^^) 

Wc{t) - —a^^\uj^)I{r,t)pb,b{t). (18e) 
Here, 5 unknown functions are subject to the following initial conditions: 

Pa,a(-00) = 1, P6,6(-00) = 7(-Oo) = /3(-CX)) = ^^^(-Oo) = 0. (19) 

It should be pointed out that only real quantities enter equations (18a)-(18e). This 
circumstance may be advantageous, especially for numerical calculations, since most 
standard numerical packages are not applicable direct to systems of ODE involving 
complex-valued coefficients. Apart from that, (16a)-(16c) and (18a)-(18e) generate 
identical results, as has been proved by extensive numerical tests carried out by us for 
laser signals with various forms of the chirp, pulse envelope, and spatial distribution. 
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However, before we demonstrate appropriate results (see section 3), it is instructive to 
present here an approximate but reasonably accurate analytic solution for C{x). This 
consideration seems to be rather useful for a qualitative description of the process as 
a whole, since it reveals some key features which make their appearance due to the 
non-stationarity of both the frequency and the amplitude of the laser signal. 

2.2. Some analytic consideration: equation (16a) 

Let us assume, in agreement with (3), that the intensity of the incident radiation is 
Gaussian both in time and space and use the most general form of the chirp, (pit/r) (cf. 
equation (4)). 

To start with, consider Riccati-type equation (16 a) to determine the complex- 
valued function C{x) : \C{x)\ = \cb^a{t)\/\ca,a{t)\: defined by (13). On introducing the 
dimensionless variable x = t/r, equation (16a) takes the form 



C"(x)-iA,„Xe ^C\x)-i 



4:71 AuT + 2(j){x) + lADe 



C{x) 



+ iDb,aIc -'^'i =0, (20) 
where the following notations were introduced for brevity: 



—T/(r, t) =Jexp <^ -—^ V , Z=\\— -exp 



c'^^-'^' " 1^ 2al J ' V 7rc2 (jta^ 1^ 2al J ' (21) 

AD = Dan - A,6- 



Note that 



X=^/2^^ r"/(^,t)dt = 8.805 -lO^^iV^^ exp 1-^1 a.u. 
cat i_oo [ajlmmy { 2af } 

depends on the time-independent parameters of the pulse only, but not on x. 

Due to its non-linearity, equation (20) admits only numerically solving, except for 
the trivial case, Au = (t){x) = 0, where C{x) is of the form 

CM ^ ^{^i5tT^ ta.h[-i:^yi5j:T^ (l + erf (^)) 

= |1-AL>|. (22) 
T > 




However, the problem can be simplified considerably by linearizing equation (20). This 
approximation, which the consideration developed below is stemmed on, is justifiable if 
|C(a;)p ^ I (^(2^) I ^ 1 holds uniformly for —00 < x < +00 and arbitrary Au, 5 and r. 
To satisfy the latter condition, it is sufficient (but not necessary) to assume that 

X=\ID,,aM\<l, (23) 

as readily follows from equation (25). If |a) = Is, \b) = 2s and r = 0, for example, 
then A ~ A^\D2s,is{^l)\/o''^ ^ 1- This imposes the following restriction on the pulse 
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A„=4mJ (X=0.65) ,^ 

: A^=10niJ (X=1.63) / ', 

: -— A„=20mJ M.26) • 

— A,=40mJM.5) ! 




Figure 2. \C{x,Ai/,S,Aui)\ as a function of x, for various pulse energies, A^, 

and fixed values of the chirp, S = 500 Mrad • s"-'^, and the frequency detuning, 
Ap = — (5/(47r) = — 250/(27r) MHz. The laser field frequency is assumed to be of 
the form: = i/l + i'chirp(i) = I'L + ^ (1 + erf(t/T)); the rest laser parameters used 
are: at = 51 ns, r = 120 ns, r = 0. 



energy: < 6 mJ, where particular values of the atomic matrix element, £>2s,is(<^l) 
(see table 1), and the typical laser parameters have been used. The above estimates are 
illustrated in figure 2 where |C(a;, Az/ = —S/4:']T,S,Ai^)\ obtained by direct numerically 
solving equation (20) is plotted versus the dimensionless variable x, at several A^-values. 
In addition, for each pulse energy, the appropriate A-value defined by equation (23) is 
indicated for reference as well. 

The linearized form of equation (20) reads: 



C'{x) - i 



AttAut + 20 (x) + lADe 



C{x)+iDh,aIe ^ =0. (24) 



Its solution obeying the zero initial condition (17) at a; = — oo has the form 
C{x) = -iDb,aIex.-p jix (AtttAu + 2(i){x)/x) + i^^^XADerf (^^^^ 



/ exp <^ - 

J — oo L 



at ( TU \ r^u 



2, ,2 



X / exp <^ -iu (AnrAu + 2(l)(u)/u) - iJ- — IADeri - \ Au. (25) 



2 r Vy2cTt7 2a. 



To demonstrate the fair accuracy of this relation for A^^ = 4 mJ and 0(,t) = 
(5rx' (1 + erf(a;)) /2, for example, the graph of \C{x, Au, 6, A^)\ ~ |c2s,is(^)|/|cis,is(t)| 
as a function of x is presented, for several chirp values, in figure 3, along with 
appropriate graphs for numerical solution of the non- linear equation (20). Together 
with plots in figure 2, these graphs provide instructive information about the magnitude 
of |c2s,is(^)P/|cis,is(i)P at various time moments. In addition, both figures clearly 
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-li -1.0 -0.5 0.0 0.5 1.0 1.5 
x=t/x 

Figure 3. \C{x,Ai',6,Aui)\ and \C{x,Al',5,Aui)\ as functions of x, for fixed values of 

the pulse energy, A^^ = 4mJ, frequency detuning, Az/ = — (5/(47r) MHz, and for several 
chirp values: 6 = 0, 10, 100, 500Mrad-s~^. The frequency of the non-monochromatic 
laser field is assumed to be of the form: u{t) = ui,-\- VcMipit) = J^l + ^ (1 + erf(t/T)), 
and the rest laser parameters used are: at = 51 ns, r = 120 ns, r = 0. 

numerical solution of equation (20); analytic solution, [(^(x)!, of equation 

(24), as given by (25). 

demonstrate a rather peculiar way in which the values for the above ratio are (strongly) 
distorted due the presence of the chirp. 

It is relevant to note that C{x) approximates the function C{x) qualitatively 
correctly for even those Ai_^ which violate condition (23). An elaborate comparison 
with exact numerical solutions of (20) shows that, for most applications, the region of 
the validity of equation (25) can be safely extended up to < 20 mJ, depending on the 
rest laser parameters. At higher pulse energies, the contribution of non-linear effects 
becomes so essential that the neglect of the ~ C^(a;)-term is no longer permissible. This 
can be seen, for example, in figure 2 where, for \x\ < 1 and Ai_j — 20 mJ, appropriate 
non-linear contribution amounts to 20%. 

As follows from the method of our derivation, equation (25) can be readily 
generalized on a wide class of similar multi-photon ionization processes that are induced 
by not very intense laser fields and which involve the chirps and laser amplitudes, such 
that \(j){t)\ <^ LOi, and I{r,t) — > sufficiently fast as i — > ±oo. Note that the latter 
condition ensures the rapid convergence of the integral in (25). In addition, it makes 
certain that the values of all physical quantities of interest are determined actually by 
only a somewhat narrow domain of x within which the laser intensity is peaked. 

Equation (25) is suitable to be used further for the analytical study of the 3-photon 
ionization with muonium, where \a) = Is, \b) = 2s and (f){x) = 6tx (1 + crf(a;)) /2. In 
particular, C{x) readily enables one to obtain analytic formulae for the probabilities: 
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|cis,is(^)P, |c2s,is(^)P, and W^{Ai/, A^, r: [0]). However, these expressions are not given 
explicitly here as they turn out to be too bulky. Instead, the formula (25) will be used 
in the following subsection 2.3, while addressing an interesting complementary problem: 
a calculation of the integral of the 3-photon ionization line profile over the entire range 
of the frequency detunings, — oo < Au < +oo. 



2.3. Further analytic consideration: the integral over the frequency detunings 

In what follows, it is shown that, for given and r, the integral of the 3-photon 
ionization line profile, iyg°°(Ai/, yl^^, r; [0]), over the entire range of the frequency 
detunings Ai/, 

/+0O 
W:^iAiy,A^,r;[m^^, (26) 
-oo 

is almost independent of particular form of the chirp, at least for low pulse energies. 
To be more precise, equation (26) defines in fact a non-linear functional of 0(i) whose 
"strength" depends, though, on the rest laser parameters. Within certain range of 
these parameters {A^, in the first instance), the functional relation (26) turns out to 
be weak, so that p{A^,r; [</)]) plays in this case the role of a sort of adiabatic invariant 
for 3-photon resonant ionization processes occurring under the action of the fields with 
slow-varying frequencies. Hence, an "almost-conservation" of p(A^,r;[0]) provides a 
simple and easily verifiable approximate criterion which enables the photoionization 
line profiles corresponding to different values/forms of the chirp to be compared and 
normalized. 

The problem discussed here was first addressed, to our knowledge, in [32], within 
the framework of the similar 3-level model and for low pulse energies only. In contrast 
with our present approach, however, the genuine photoionization rate of the level \b) 

(that is, 25") was set in [32] equal, respectively, to zero and to some time- and energy- 
independent constant, while describing the dynamics of the 1st and the 2nd stages of the 
resonant 3-photon ionization. Although this approximation provides acceptable results 
for < 6mJ, say, it is of interest, both experimental and theoretical, to reconsider the 
entire problem by releasing the above simplifying assumptions and getting, thereby, a 
deeper insight into the physics of the process. 



2.3.1. Low A^^ Let us consider first the simplest case of low A^, such that condition 
(23) is satisfied. In addition, we will still assume, without a loss of generahty, that the 
intensity of the incident radiation is Gaussian both in time and space, without specifying, 
though, particular form of the chirp for a moment. By making use of equation (16c), 
this allows p{A^,r] [(j)]) to be expressed in terms of C{x) as 



p{A^,r-[<p]) 



+ 00 



dAv 



X exp I 2Db,aI^ 



/X 
oc 



+ 0O 



dxexp 



C{x,Aiy, A^,r] 



exp 



2at 



C{u, Av, A,^, r)du 



(27) 
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Here, our former notations (21) have been employed. To this end, it must be noted that 
IAtttAu + 20(0) + XAD\ > 1 is satisfied for all Ap and X, except those Au which are 
in the vicinity of the lineshape's maximum, Au ~ — </)(0)/(47rr), and for A^^ such that 
\IAD\ < 1. The latter condition happens to be more restrictive as compared to the one 
assumed here. This is due the fact that the diagonal matrix elements are usually much 
bigger than the non-diagonal ones (see table 1 for comparison). By making use of (25), 
(5) and presuming that the above condition is fulfilled, one can show that 



2D, 



b.a- 



exp 



C{u, Ai/, Ai_j, r)du 

a(AD) 



47rTAi/ + 20(O) + JK(AD) +J2[$>(AD)]2 



1 + erf 



TX 



V2i 



(28) 



Note that, for any x, arbitrary laser parameters and D^ a, this expression is negative, as 
it should be, since ^{AD) = — Q=(Db^{,) = (47ra)~V^'''''(a;L) > 0. For small pulse energies 
considered here, the right-hand side of (28) is proportional to I^, which permits the 
appropriate exponent in (27) to be expanded in terms of A^. On retaining two leading 
terms in this expansion, this yields 

'^^^ ' ' ~ d Ai/ / da; exp 



p{A^,r;[^])^^-ai'\u;^)I 



T^X^ 



2at 



C{x,Aiy, A^,r) 



X 



27rDlJ'<^iAD) 



T 



1 + erf 



TX 



\f2. 



0-t 



(29) 



47rTAi/ + 20(O)+X3?(AD) +X2[^>(AL>)]2 

Further calculations are straightforward but cumbersome. Hence, some technical details 
will be given below for the first term in the curly brackets in (29) only; the corresponding 
contribution is denoted as p'^^^^A^^ r; [0]). In this simplest case, by substituting (25) for 
C [x] and interchanging the order of integration, one arrives at 



C (7, 



47r 



exp 



t'^x^ 



2at 



X 



exp + ^/2^^TS(Ar>)erf i 



TU 



VV2(7t) 



dw 



47r ^T^ui^ %AD) { 



exp 



-2u' - V27r—IQ{AD)eric{u] 



T 



du 



(30) 



For two most important cases of interest, \I'<s{AD) \ <^ 1 and \I'^{AD)\ > 1, the integral 
in (30) can be evaluated, respectively, by expanding the integrand into the power series 
in I'^{AD) and by using the method of the steepest descent. This yields 

+0O 



exp 



-2ii^ - V2n—IQ{AD)eric{u) 



du 



- 7i{at/rMAD)I 
(V^rV23/2a2) mAD)T 



-2 



for \I'^{AD)\ < 1 
for \m{AD)\ > 1 



(31) 



Resonant 3-photon ionization 



15 



Accordingly, we are left finally with 

^n2 ^/ A/^(^t/r)X=^ for |X$>(Ai?)| « 1 

- 47r V 2 t^cl'"' ^"""^ \ [v/2S(AD)]-iX^ for |J9(AD)| » 1 ' ^'^ ^ 

Note that in either limiting case, the result is independent, in accord with our initial 
statement, of particular form of the chirp 0. To clarify this point, one should emphasize 
that the only two simplifications made in the course of derivation were: (i) the use of 
C{x) given by (25) and (ii) the replacement of the appropriate exponent by a unity. Both 
approximations are fully justifiable for : \A^Di,^a{^L)\ ^ 1 assumed here. Obviously, 
the latter condition is consistent with the first relation in (32), which has been formerly 
obtained by means of a different technique (and cast in a slightly different analytic 
form) in [32]. Apart from being chirp-independent, this result demonstrates also that 
p^'^\A^,r;[(l)]) oc A^, as opposed to the next limiting case, that is, \I'^{AD)\ > 1, 
where p*^°^(At^,r; [(f)]) oc A^. A gradual decrease of the exponent is clearly evidenced 
by numerical simulations carried out by us for muonium, while scanning sufficiently 
broad A^-domain. Some results of this numerical study are discussed in more detail 
in section 3 (see figure 8, for example). These show, in particular, that equations (32) 
describe p^^^A^^^r; [0]) quahtatively correctly for both low and moderately high pulse 
energies. 

In a similar manner one can also obtain appropriate contribution, denoted 
p^^\Ai^,r ;[(()]), coming from the second term in the curly brackets in (29). After 
tedious calculations whose details will be given elsewhere, one ends up with the following 
expression for a sum of two contributions: 

pW(A^,r;M)+p«(A^,r;M) 

A"^ M^la^ (l - ^^S^M^) for \m{AD) I « 1 



„2 

C ~ 

C VTT CTt (7j 

47r 2 T^^b 



M^I" (l - '-^$011^) for \m(AD)\ > 1 



Here, according to (5), we have discarded the terms in p^^\Ai_j,r; [0]) containing the 
second and higher order derivatives of (p{x). This simphfication results in the fact that 
entire equation (33) turns out to be completely independent of the chirp. 

The right-hand side of equation (5) depends (through X) on the radius r. To 
eliminate this dependence and to facilitate, thereby, an adequate comparison with 
the results of numerical calculations, the left- and right-hand sides of (5) can be 
integrated over the entire XOF-plane perpendicular to the direction of the laser beam's 
propagation (2;-axis). A clear physical meaning of such an averaging procedure is 
discussed in the comments following equation (36). 

By virtue of the identity, 



/ + CO P 
■oo J — 
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Figure 4. The integral of the spatially-averaged two-step 3-photon ionization 
probabilities of the ground state of muonium, {W^ (Ai/, 6, A^,)) , over the entire range of 

the frequency detunings Av, as a function of the chirp's magnitude, S: three particular 
pulse energies are considered: A,^ = 4, 10, 20 mJ. The chirped frequency of the laser 
field employed is of the form: ^{t) = ul + i'chirp(i) = + ^ (1 + erf(i/T)). The rest 
laser parameters used are: = 51 ns, r = 120 ns, (Tr = 0.64 mm. 



which defines the spatially-independent part of X, denoted as X = •\/2/ nc^TA^/ ata^ , 
this yields eventually the following r-independent result: 

m « / / {p^'\A^,T- [0]) +p«(A,,r; [0])) dxdy 



t^"^-^\^^)^la^ (l - '-^'iKa^) for \X^{tAD)\ « 1 

f34) 

t'^'^y^\^^)^?'{}-'^iKa^) forix$>(Az^)i>i ■ 



In order to check this equation, the numerical calculation of p{^A^^\ [0]), in a wide 
range of and for typical values of the rest physical parameters used in the 1S-2S 
experiment in muonium, have been carried out. This was done by an extra numerical 
integration of the 3-photon ionization line profiles (W^c°^(Az/, 5, A^^)) (see equation (36) 
for definition) over the entire Az/-domain. These results are shown in figure 4, for 
three particular pulse energies. At ^4^^ = 4 mJ, being the only value displayed 
which satisfies condition (23), the plot demonstrates that p(At^; [0]) is conserved to 
the 1% accuracy, within a wide range of the chirp's magnitude. The mean numerical 
value, p{A^ — 4 mJ,(7i. = 0.64 mm; [0]) 0.080 MHz, is to be compared with that 
obtained by means of the f\,rst equation in (34): p(Aj^ = 4 mJ, (Tr = 0.64 mm; [0]) 
0.1242 — 0.0292 = 0.095 MHz. To make these results comparable, the latter value 
includes also an additional weighting factor, (97r(Tr)~^. Even though parameter A is 
very close to unity: A = \X'^{D2s,2s)\ — |^-D2s,is| ~ 1, for A^^ = 4 mJ and 0"r = 0.64 mm 
used here, numerical and analytical data are consistent to 20% of relative accuracy. 
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Interesting to note that the first /second relation in (34) overestimates/underestimate 
the true value of p(A^; [0]). Much better agreement can be anticipated (and actually 
takes place) for lower pulse energies and/or bigger spatial dispersions a^. For = 2mJ 
and CTj. = 1.5 mm (apppropriate line profiles are displayed in figure 5(a)), for example, 
numerical and analytical results read: p{A^ = 2mJ, = 1.5mm; [0]) pa 0.90 ■ 10~'^MHz 
and (0.936 - 0.02) • 10"^ = 0.916 • 10"^ MHz, respectively. 

2.3.2. High For higher pulse energies and/or smaller spatial dispersions (7r, such 
that \iQ{AD)\ » 1, equation (34) fails to describe adequately the true function 
p{A^^; [(f)]). This refers both to the absolute values of p{A^; [0]) and, particularly, to 
its dependence on the magnitude of the chirp, which makes its appearance at high 
pulse energies. It has been elucidated above that the main underlying reason for this 
failure originates from the fact that "the linear approximation" ceases to be valid any 
longer, due to corrections caused by the C^-term in equation (20). This non-linear 
term does not allow any reasonably accurate analytic treatment to be developed, so 
that the numerical analysis of the problem must be used instead. To simulate the 
high-energy regime, the results of our calculation are presented in figure 4, for two 
pulse energies, Ai_^ — 10, 20 mJ, and for typical values of the rest physical parameters 
used in the 1S-2S experiment in muonium. As in the case of low A^^, this was done 
by the numerical integration of the 3-photon ionization line profiles {W^{Au,S, A^^)) 
over the entire Az/-domain. The plots demonstrate in fact a noticeable deviation of 
p(A^; [0]) from a constant value, with an increase of the pulse energy. For example, 
at = 20 mJ and for 6 in the range 6 = 0... 500/ (27r) MHz, this deviation (that is, 
the dependence on the chirp's magnitude) amounts already to 30%; in absolute units, 
the data read: p{A^ = 20 mJ, (7^ = 0.64 mm; [0]) = 2.265 MHz for (5 = MHz and 
p(A^ = 20 mJ, (7r = 0.64 mm; [0]) = 2.90 MHz for 5 = 500/(27r) MHz. Furthermore, the 
gauge of the given deviation may be expected to be even more pronounced at higher pulse 
energies, thus preventing us actually from viewing p{Ai^; [0]) as a conserving quantity 
in this case. Although such A^^ are of no practical interest at the moment for the lS'-25' 
experimental study with muonium, the entire effect might be of relevance for future 
experimental studies where much higher laser intensities are involved. 

3. Application to the experiment in muonium 

Let us apply the results obtained in preceding sections to particular laser parameters 
adopted in the 1S-2S experiment in muonium [8], by assuming that \a) = IS, \b) = 2S 
and setting (pit) = St (1 + erf(t/r)) /2. It is convenient in this section to measure the 
energies in the units: h = = m* = 1, where m* = md (1 + mc/m^) ^ (1.00484)^^ mg 
stands for the reduced mass of the electron and muon. Here, the value for the electron 
to ^~ mass ratio, m^/m^ = 206.768262, was used. The quantity a* = t^/{m*e^) = 
{me/m*)ao stands for the unity of distance, with = /{irif.e'^) = 0.529 • 10~^ cm 
being the Bohr radius. Accordingly, one-particle energies and dynamic polarizabilities. 
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Table 1. The values of the two-photon matrix elements -Db.a, (b,a) = Is, 2s, defined 

by the sums in equations (10)-(12), and the value of the 1 photon photoionization 
cross section of the 25* level, given by (35), at cjl ~ 3/16 a. u. (Al = 244 nm); linear 
polarization of the photon, e, is assumed in both cases. 

Dls,ls{i^h) D2s,2sii^h) -D2s,1s(wl) ^^^^(wl) 

-r,.7141 29. 8r,.3r, - i 12.8232 7.8.^).3.^) 0.220-5 



a^^{uj), are measured in the units of e^/a* and (a*)^. 

As was mentioned above, 3-photon ionization probabihty, W^{Ai', S,A^,r), can be 
obtained by accurate numerically solving either the system (166)-(16c) or (18a)-(18e). In 
both cases this has been done by the stepwise time integration, starting from to — —Sat, 
while scanning sufficiently broad Ai/-domain centered at Ai/ — (i.e. cui, — ^0U2s,is), 
for all r = \/ x'^ + y'^ on the r-mesh: = ka^-flO, k = . . . 50. Numerical values 
of the matrix elements, Dh a, {o-, b) = Is, 2s, which enter (166)-(16c) were calculated 
in [30, 17, 31]; these arc compiled in table 1, together with the value of the single 
photoionization cross section of the 25'-level at a; = a;L. The latter is given explicitly 
by (see [31] and references therein) 



Ol4_2 / T \ / T \ * „-4T?arctan(2/n) I A j 



^'^s'M - ^« 1 + 3^ pi rj = \ (35) 

where a is the fine structure constant and = 1/8 denotes the ionization potential of 
the 2S'-state. 

At the final stage of our numerical procedure, a spatial averaging has been 
performed in order to obtain r- independent ionization profiles, {W^{Au,S,Ai^)). This 
has been achieved by an extra integration of W^{Ai', 5, A^, r) over entire XOF-plane 
perpendicular to the direction of the laser beam's propagation: 

kWT{Av, S, A^)) = ^ y J W^{Au, S, A^, r) dxdy. (36) 

Basically, this newly introduced quantity can be interpreted as an averaged ionization 
probability related to the entire laser beam spot, since (W^lAu, 6, A^^)) is independent 
of the distance from the beam's axis. Alternatively, in view of the low values of the muon 
density in the media [8] (see below), {W^{Ai',d,A^)) can be called as "the ionization 
probability per one (muon) atom in the beam" . 

An auxiliary factor, l/{97ia^), has been introduced in (36) in order to make 
{W^{Au,6,A^)) dimensionless, as the original probability W^{Ai/,5,A^,r) is. This 
particular choice was adopted according to the so-called "3(T"-rule being inherent to 
various problems involving Gaussian law. In the case considered, the "3(7" -rule makes its 
appearance through the fact that the values of W^{Ai',6, Ai_j,r) happen to be almost 
negligible for r > 3crr, unless the pulse energies higher than A^ = 20 mJ, say, are 
considered. For such a high energy regime, the probability W^lAu, 6, A^, r) has a long- 
distance "tail" spreading out beyond the effective beam radius, r = 3crr, within which 
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the values of W^°°(Az/, S, A^, r) may be very close to, yet less than, 1. In this situation 
the proper normalization factor different from that chosen above must be used so as 
to prevent (lV^°^(Az/, 5, A^^)) from being bigger than 1. As is clearly evidenced by our 
calculations, the values of (lV^°^(Az/, 5, A^)) defined by equation (36) exceed 1 starting 
already from A^ ~ 20 -j- 25 mJ, depending on the rest laser parameters. One should 
point out, however, that such are hardly attainable currently with proper pulsed 
lasers sources in the required frequency range. This makes the high energy regime in 
the resonant 3-photon ionization with muonium to be mostly of academic rather than 
practical interest at the moment, thus justifying the definition of equation (36). To 
avoid any misunderstanding, it must be clearly stated once again that the spurious 
effect mentioned here originates solely from the particular choice of the normalization 
factor in (36), and has, therefore, nothing to do with either the model developed in 
the work or a lack of the numerical accuracy employed by us. The latter was chosen 
to be equal to 10~^ which proves to be sufficient to ensure that W^{Ai',d, A^^,r) < 1 
holds (as it should!) for any particular individual set of all parameters involved in the 
problem, including r, Au, 5, r, and (7r. 

The averaged probability {W^{Au,S, A^^)) readily enables one to estimate the 
expected number of experimentally observed ionization events, A^event) defined as the 
number of the muonium atoms which are ionized by a sequence of Aghot identical 
Gaussian laser pulses (3) within the time T and detected eventually in apparatus. Such 
a formulation of the problem corresponds direct to the actual experimental situation 
where the gas media containing muonium atoms interacts with Ashot = 25 laser pulses 
per second. Under these conditions, the required number of the ionization events to 
be detected during the observation time T » AghotT" can be estimated by means of the 
following simple relation: 

Here, po — 4 ■ 10^'^ atoms/mm^ is the spatial density of the muonium atoms in the 
chamber, S = Qvrcr^ ~ 11 mm^ and L ^ 10"^ mm stand, respectively, for the effective 
cross section of the laser beam (see above) and the total path which this beam travels 
in the media, rj ~ 10% denotes an efficiency of registration of the muonium atoms 
in apparatus; finally, W^^{S,Ai^) is the value of the averaged 3-photon ionization 
probability {W^{Al',5,A^^)) at its maximum, for particular values of the chirp and 
the energy of the pulse. Provided S and are fixed, this maximum is achieved 
approximately at Au ^ —6 /{An), almost independently of A^j, as can be seen in 
figure 5(b). In particular, W^'''^{5 = 10 MHz, = 4 mJ) ^ 0.014 thus leading, for 
example, to the following expected number of events detectable during the observation 
time T = 3600 s: Aevent — 5800. It is rather instructive that this value happens to 
be quite close to that found preliminary in the IS — 2S experiment in muonium being 
currently underway. 

Some results of our simulations are presented in figures 5-9, for various magnitudes 
of the chirp 5 and various parameters of the laser pulse which are typically used in the 
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measurements of the IS — 2S energy separation in muonium; particular values employed 
are indicated in figures 5,6. Both these graphs demonstrate appreciable red-shift of the 
maxima (i.e. towards lower frequencies, relative to z/l = u;2s,is/47r) of the photoionization 
profiles with an increase of the chirp. It must be noted, however, that this shift is 
caused by a combination of two competing mechanisms: (i) the shift being due to the 
combined AC Stark effect for the IS- and 25'-levels and (ii) that arising because of the 
chirp itself. The curves show that the former mechanism turns out to be small compared 
with the latter at relatively small pulse energies, A^^ — 2 . . A mJ. At A^^ — A mJ and 
S — 0, for instance, i.e. in the case of completely unchirped laser signal, the AC Stark 
shift amounts approximately to 1 MHz. One can anticipate, however, that the relative 
contribution of the AC Stark shift will increase as the laser intensity increases, since 
its magnitude is linearly proportional to the laser intensity, whereas the chirp-induced 
shift is almost intensity-independent. This behaviour can be seen in figures 6(a,b) and 
is summarized in figure 7(a) showing, in particular, that the AC Stark shift becomes 
equals to 4 MHz at A = 20 mJ, S — 0. In addition to the shift of the maxima, the 
curves in figures 5,6 exhibit a strong dependence on both the spatial dispersion a^. and 
the power A^ of the laser signal. It should be noted that an appreciable distortion 
of the photoionization line shapes happens to be particularly enhanced at relatively 
large chirp values, where the line profiles become asymmetric. Generally speaking, this 
asymmetry is an intrinsic feature of the 3-photon resonant photoionization occurring 
under the action of the chirped pulsed laser signal. Moreover, an asymmetry of the 
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Figure 5. Two-step 3-photon ionization probabilities, {W^{Ai',6,Aui)), of the 
ground state of muonium by the laser field with the chirped frequency ^{t) = 
i^L + J^chirp(i) EE + ^ (1 + erf (t/r)) , S = 500, 300, 100, 70, OMrad • s-\ versus laser 
frequency detuning, Av = (wl — ^0J2s,is) /27r. Two pairs of different values for the 
pulse energy, A^^, and the spatial dispersion of the laser signal, a^., are used: (a) 
A = 2 mJ, (Tr = 1.5 mm, and (b) ^4 = 4 mJ, ar = 0.64 mm. 
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Figure 6. Two-step 3-photon ionization probabilities, {W^{Ai/, (5, A^)), of the ground 
state of muonium by the laser field with the chirped frequency = 1/1, + fchirp(i) = 

t'L -I- ^ (1 + erf(t/T)) ,5 = 0, 100, 300, SOOMrad • s~^, versus laser frequency detuning, 
Av = (wl — ^C02s,is) /27r, for two pulse energies: (a) A,^ = lOmJ and (b) A,^ = 20mJ. 



ionization line shapes happens to be one of the most pronounced manifestations of the 
non-monochromaticity of the laser field. This phenomenon is basically due to the fact 
that the chirp, even though its relative magnitude amounts to only ~ 10^*^ . . . 10^^ 
in our calculations, violates an equivalence of initial {t — —00) and final {t = +00) time 




Figure 7. (a) - the frequency red-shift of the photoionization profile's maximum, 

against the pulse power A^^. The chirped frequency of the laser field is assumed to be 
of the form: i^{t) = i/L+i^chirp(t) = 1^+^ (1 + erf(Vr)) ,6 = 0, 100, 300, 500Mrad-s-i; 
(b) - the power-dependence of the line profile's maximum value. In both cases the rest 
laser parameters are the same as in figures 6(a,b). 
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Figure 8. The A^-dependence of the exponent, /x(A^j), defined by the relation: 



moments. Indeed, the photoionization probability turns out to be either more efficient 
or suppressed, depending on whether the time-dependent chirped laser frequency is in 
the resonance with the 2-photon \a) — > \b) transition or slightly off it. The process 
as a whole becomes then somewhat "time-sensitive", although the form of the laser 
pulse envelope used is time-invariant. It should be noted, however, that the gauge of 
the above asymmetry depends strongly on the laser power, as well as particular forms 
of the laser amplitude and the chirp. This dependence can be seen in figures 6(a,b) 
where higher A^-values as compared with figures 5(a,b) are used. The former show 
that the asymmetry discussed is less pronounced as long as higher laser intensities are 
involved. Two essential characteristics of the line profiles' maximum, that are, the 
shift and the magnitude, are plotted in figure 7(a,b) as functions of the laser power for 
various values of the chirp. In particular, the graph 7(b) provides information about 
the dependence of the fine shapes' maximum value, W^^{5,Ai^), on the pulse energy 
A^^. It is quite natural to parameterize this relation by the following simple power law: 
W^'^{S,A^) = 7((5)A'^(^-). The exponent here proves to be a slow-varying function of 
the pulse energy, such that {dfi{A^) / dA,^) / ii{Ai^) ^ (Ai^lnAi^)'^. According to this 
condition, /^(^a,) can be estimated as 



One should note, for example, an appreciable deviation of the fi{Ai^)-values from 3 
which must be expected for a 2-step 3-photon ionization of an atom in the state \a) by 
a weak, monochromatic and spatially homogeneous laser field, without the relaxation of 
the intermediate levels. For 6 — 0, the A^^-dependence of IJ,{A^) is displayed in figure 8. 
This graph demonstrates that Ijl{A^ — 2 mJ) Ri 2.75 gradually decreases as the pulse 
energy increases and tends, for ^ 20 mJ, to an almost constant value, /i ~ 0.5, thus 
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Figure 9. Time-dependence of the chirp - (a) - and the laser power - (b) - 
both detected in the former experimental measurements of the 15 — 25 separation in 
muonium [8]; (c) - photoionization line profile measured in [8], versus laser frequency 
detuning, Az/ = (wl - ^i^2s,u) /^tt- 



indicating on the presence of saturation in the 3-photon transition considered. 

Finally, it should be mentioned that the above consideration has been developed 
mainly in the attempt to simulate those line profiles of the two-step 3-photon ionization 
probabilities that are supposed to be measured soon in the new IS — 2S experimental 
investigation of muonium. Essential physical parameters employed in this study are still 
not fixed completely and are subject to further changes. It is therefore rather tempting to 
apply our present technique to appropriate data formerly collected within the framework 
of the previously employed experimental setup [8]. These are shown in figures 9(a,b), 
along with the photoionization line profile obtained within the framework of the current 
theoretical approach, figure 9(c). The fact that the corresponding experimental line 
shape (not shown) is virtually indistinguishable from the theoretical one provides an 
additional and independent check of the overall validity of the model developed in this 
work. 

4. Conclusion 

We have developed a simple theoretical scheme intended to describe, to the 1 MHz 

accuracy, the stepwise 3-photon resonant photoionization in hydrogenic systems, 
induced by the chirped laser field with time-dependent amplitude. It has been shown 
that such an accuracy can be achieved within the framework of a simple 3-level model, 
by taking into account (i) the AC Stark shifts and (ii) non-zero ionization rates of the 
levels involved, together with (iii) a spatial inhomogeneity of the laser signal and (iv) 
arbitrary ^-dependencies of its intensity (pulse envelope), I{r,t), and frequency, Lj{t). 
The system of equations (16&)-(16c) or, equivalently, equations (18a)-(18e) is a key point 
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of the method employed. These are of independent significance as being not specifically 
restricted to particular states of reference, laser intensities, and chirps, so that the results 
obtained for IS*—, 25*— and eP-states can be generalized on arbitrary \a), \b), \c) levels of 
hydrogen-like atoms. Excited ns— levels, such that 2 < n < 5, are of particular current 
interest for ultra-high precision laser spectroscopy, and it should be expected that an 
adequate interpretation of experimental data should require an accurate theoretical 
account of a wide spectrum of light-induced effects. Some of these effects are beyond 
the scope of our present consideration and, in the first instance, comprise appropriate 
corrections due to inevitable motion of atoms in the media and associated second order 
Doppler shifts, as being of major importance. Indeed, even for CW lasers where the 
signals are usually almost unchirped in the laboratory frame, the amplitude and the 
frequency of the pulsed laser signals become essentially time-dependent in the atomic 
rest frame, thus leading, as was demonstrated above, to the shifts and distortions of 
ionization/excitation lines. The relative contribution of the effects discussed is estimated 
as a'^v^ with Vz being the atomic velocity in the laboratory frame, and these should be 
incorporated in the theoretical scheme developed here, along with relativistic/radiative 
corrections ~ to the operator of the particle-laser field interaction. This work is 
currently in progress, and we consider it as a subject of forthcoming publications. 
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